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ABSTRACT 

We report results of RXTE observations of PSR B1951+32 using the 
PCA instrument for 19k seconds during 1996 March 24th. We apphed the 
contemporaneous radio ephemeris and various statistical tests to search for 
evidence of pulsation. These analyses yield intriguing yet inconclusive evidence 
for the presence of the pulsation in the time series: confidence level for the 
presence of pulsation is 94.5% in the 2.0-4.8 keV band and 97.6% in the 
4.8-6.3 keV band based on the H-test. Under the premise of non-detection 
of pulsation, we derive estimated 2-a upper limits for the pulsed flux to be 
3.3 X lO-^cm-^s-^keV"^ in the 2.0-4.8 keV band, 2.7 x lO-^cm-^g-ikeV"^ in 
the 4.8-8.5 keV band, and 2.0 x lO-^cm-^s'^keV"^ in the 8.5-13.0 keV band. 
These upper limits are consistent with the trend of spectral turn-over from 
high-energy gamma-ray emission as suggested by the OSSE upper limit. Such 
turn-over strongly suggests the outer magnetosphere as the emission site for 
pulsed gamma-rays. These RXTE upper limits for X-ray pulsation are, on the 
other hand, not consistent with the extrapolation of reported power-law spectra 
from the point source observed by ROSAT in the 0.1-2.4 keV band, assuming 
a constant pulse fraction: The pulsed soft X-ray emission detected by ROSAT 
must follow a much softer spectrum than that of the overall point source. 

Subject headings: pulsars: individual: PSR B1951+32 
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1. Introduction 

The pulsar PSR B1951+32, located at the center of the morphologically peculiar radio 
nebula CTB 80, is a 39.5-msec radio pulsar (Clifton et al. 1987; Kulkarni et al. 1988) with 
a characteristic age of 1.1 x 10^ yr and an inferred surface dipole magnetic field of 4.9 x 10^^ 
G. An X-ray point source was observed within the X-ray nebula related to the radio core 
of CTB 80 (Becker, Helfand & Szymkowiak 1982; Seward 1983; Wang & Seward 1984). 
Search for X-ray pulsation from this source with EXOSAT yielded inconclusive evidence for 
pulsation (confidence level of 97%, by Ogelman & Buccheri 1987, and 93% by Angelini et al. 
1988). The pulsed emission was detected by ROSAT at a 99% confidence level (Safi-Harb, 
Ogelman & Finley 1995), which shows a single peak roughly consistent in phase with the 
radio emission. The overall spectrum of this point source in 0.1-2.4 keV is best fitted by a 
power law with a photon spectral index —1.6 ± 0.2, and an estimated pulsed fraction of 
0.35. 

The EGRET instrument on CGRO observed and detected gamma-ray pulsation from 
PSR B1951+32 above 100 MeV (Ramanamurthy et al. 1995), making it a member of 
EGRET hard gamma-ray pulsar family with a similar age to Vela and Geminga. The 
gamma-ray lightcurve shows two peaks at phase 0.16 and 0.60 with phase 0.0 being the 
radio peak. Its spectrum in the EGRET energy range follows a power law with a photon 
spectral index of about —1.8 (Ramanamurthy et al. 1995; Fierro 1995) over about two 
decades of photon energy. Recently, pulsed emission is reported from the COMPTEL 
instrument in the 0.75-10.0 MeV band (Kuiper et al. 1996). The OSSE and BATSE 
instruments on CGRO only reported upper limits of pulsed emission in the lower energy 
band (Schroeder et al. 1995; Wilson et al. 1992). 

There have been a number of models proposed to explain the gamma-ray emission with 
dramatically different emission sites, some at or very near the surface of the neutron star 
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and some very far away. Recently, Ho & Chang (1996) proposed a geometry-independent 
argument to constrain the possible site and mechanism of the gamma-ray emission based on 
the commonality of power-law emission in the EGRET (and possibly COMPTEL) pulsars. 
In such arguments, it is important to know whether and how the gamma-ray power-law 
spectra turn over towards low energy. (See Section 4 for more discussions.) 

To gain better understanding of the overall spectral behavior, especially between keV 
and MeV, we conducted an observation of PSR B1951-I-32 using both PCA and HEXTE 
on board RXTE during Cycle 1. Analysis of the 19k-second PCA data does not yield 
conclusive evidence for pulsation from 2.0 to 13.0 keV. The derived 2-0" upper limits provide 
support for the hard turn-over for the high-energy gamma-ray emission. It also indicates 
that the soft X-ray pulsation observed by ROSAT has a very soft spectrum. We described 
the observation in Section 2. The analyses and results for the PCA data are discussed in 
Section 3. We discuss the theoretical implications of this observation and future work in 
Section 4. 

2. Observation 

The PCA and HEXTE on board RXTE were pointed at PSR B1951+32 on March 24, 
1996 (MJD 50166), for about 10.5 hours including earth occultations. The RXTE mission, 
spacecraft and instrument capabilities are described in Swank et al. (1995), Giles et al. 
(1995) and Zhang et al. (1993) The PCA consists of five essentially identical PCUs with a 
total effective area of 6729 cm^, with no imaging capability. The field of view is one degree. 

After examining the data, two exclusions were applied to the data set. First, data from 
the PCA pulse-height channel 36-255 (13.0-90.0 keV) are excluded due to high instrumental 
noise. Second, we observed unexplained anomalous increase during two intervals of our 
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exposure. Under the advice of RXTE Guest Observer Facility experts, data obtained during 
these two intervals were excluded. In the second half of the observation, two of the five 
PCUs were turned off. The overall usable data used for this analysis contain two segments 
of 8928 s X 6729 cm^ and 10304 s x 4180 cm^ for a total of 1.03 x lO^scni^, or equivalently, 
a total integration time of 19232 seconds and an average effective area of 5363.3 cm^. 

Around the same epoch of the RXTE observation, PSR B1951+32 was also monitored 
at Jodrell Bank Radio Observatory. The radio ephemeris is summarized in Table 1 and 
used as the input for pulsation search. 

EDITOR: PLACE TABLE I HERE. 



3. Analysis and Results 

The data were reduced to the barycenter and analyzed using the JPL DE200 ephemeris, 
the pulsar position listed in Table [^, and standard RXTE reduction package (FTOOLS 
v.3.5.2 and XANADU/XRONOS v. 4. 02). Lightcurve folding was performed separately for 
each of four typical energy bands and various combinations using the radio ephemeris in 
Table |l]. The four typical energy bands are commonly designated as band 1 through 4, with 
each covering PGA channels 0-13 (2.0-4.8 keV), 14-17 (4.8-6.3 keV), 18-23 (6.3-8.5 keV), 
and 24-35 (8.5-13.0 keV), respectively. None of the folded lightcurves showed significant 
deviation from a model steady distribution under the Pearson's x^-test (Leahy et al. 
1983a,b). Specifically, the x^ values for the folded lightcurves shown in Figure [I| are, for 
19 degrees of freedom, 27.4 for band 1, 21.1 for band 2, and 8.38 for the combined band 
of 3 and 4. In addition to instrumental and cosmic X-ray background, the DG component 
is mostly likely the contribution from the ROSAT point source and its associated X-ray 
nebula. 
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EDITOR: PLACE FIGURE |1| HERE. 

To further ascertain the absence of pulsation, we performed the bin-independent 
parameter- free H-test (De Jager, Swanepoel & Raubenheimer 1989). In this analysis, all 
detected raw photons with the corrected arrival time are used. The H-test was applied to 
the data in different energy bands and various combinations. The results of the H-test all 
show high probability that the data are consistent with steady source except for band 1 
and 2. The if-values are 7.286 and 9.334 for bands 1 and 2 (both at 171=1). This yields a 
5.5% and 2.4% probability of being consistent with a steady source. Applying the straight 
Zf-test (the Rayleigh test, which is more appropriate if the underlying pulse profile is 
sinusoidal), the probability of the data being consistent with a steady source is 2.9% and 
0.9% for bands 1 and 2. 

Based on these analyses, we do not consider the null probability, although intriguingly 
small using the H-test, provides conclusive evidence of pulsation from PSR B1951+32 in 
the XTE/PCA energy band. 

The upper limit of pulsed flux is estimated following the prescription given by Ulmer et 
al. (1991). Assuming a duty cycle of 0.5 and combining bands 2 and 3 to yield a comparable 
total number of counts to those in bands 1 and 4 individually, we obtain the following 2-a 
upper limits, which are also shown in Figure |^: 3.3 x 10~^cm~^s~'^keV~^ for 2.0-4.8 keV, 
2.7 X 10^^ cm-^s-ikeV"^ for 4.8-8.5 keV, and 2.0 x 10^^ cm^^s-ikeV"^ for 8.5-13.0 keV. 

EDITOR: PLACE FIGURE HERE. 
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4. Discussion 

Gamma-ray pulsar: The current XTE/PCA upper limits provide support to the 
combined CGRO observation that the gamma-ray pulsed emission from PSR B1951+32 
follows a power law with a significant break towards low energy (EGRET and COMPTEL 
detection along with OSSE and BATSE upper limits), as indicated in Figure |[ Such 
spectral behavior is seen in the Vela and Geminga pulsars. For PSR B1951+32, the break 
energy (photon energy at which there is a significant break in photon spectral indices to, 
say, harder than —1) is estimated to be between 70 keV and 3 MeV. As noted in Ho & 
Chang (1996), this common trait could play an important role in the theoretical modeling 
of this family of "EGRET pulsars." The power law of these pulsars typically covers two 
orders of magnitude in the EGRET band with best-fit photon spectral indices in the range 
of —1.4 to —1.8. These photon spectral indices cannot be produced by a mono-energetic 
relativistic electron distribution under currently proposed radiation mechanisms. The most 
likely origin of this power-law behavior is the cooling (energy loss) through the dominant 
radiation mechanism responsible for the gamma-rays. For example, the most simplistic 
cooling model for a steady state electron distribution will yield a photon spectral index of 
—3/2 for synchrotron radiation and —5/3 for curvature radiation. The cooling spectrum 
will continue towards low energy until the electron distrbution (in energy space) is no 
longer affected by cooling: i.e. the cooling spectrum turns hard at the break energy which 
corresponds to the location and electron energy where the radiative cooling time scale is 
comparable to the dynamical time scale of the relativistic motion of the electrons. Such 
an argument allows us to constrain the radiation mechanism and, more importantly, the 
emission site, which to date remains unsettled with great bifurcation among gamma-ray 
pulsar models. Following this argument and examining various radiation mechanisms, we 
find that, for PSR B1951+32, the gamma-ray pulsations are most likely generated by 
synchrotron radiation with a typical pitch angle of 0.1 to 0.001 and emission site is about 



10^ to 10^ cm from the star: i.e. in the outer magnetosphere. 

X-ray pulsar: Safi-Harb, Ogelman & Finley (1995) reported ROSAT observation of 
PSR 31951+32 with an estimated pulsed fraction of 0.35 and an overall spectrum following 
a power law with a photon spectral index of —1.6 ± 0.2. To date, there is no published 
spectrum for the pulsed component. Safi-Harb et al. (1995) estimated the soft X-ray 
pulsation duty cycle to be 0.1. Such a duty cycle will lower the XTE upper limit estimated 
above. It is clear from Figure |^ that the XTE/PCA observation necessitates a steep 
drop-off at around 2 keV for the pulsed component. This is consistent with Safi-Harb 
et al.'s report of no pulsation near 2 keV. It is almost certain that, for PSR B1951+32, 
the ROSAT observed pulsation below 1.3 keV is separate from the EGRET/COMPTEL 
observed gamma-ray pulsation. More than likely, they are from different origin and emitted 
at different location. 

In summary, XTE/PCA observation over 19k seconds shows no definitive evidence of 
pulsation. The upper limits can be used to help constrain theoretical models. Our analysis 
does show tantalizing hints. In addition to the small null probability from the if -test, the 
peaks for bands 1 and 2 in Figure |I|, taken at face value, are separated by 0.45 in phase, 
reminiscent to that of the Geminga pulsar in 0.07-1.5 keV (Halpern & Ruderman 1993). A 
long exposure on XTE, e.g. 100 ksecs, will provide better statistics and help advance our 
understanding of this and other similar pulsars. 

We thank the XTE team in XTE/GOF at GSFC, especially James Lochner, for their 
help in data reduction and analysis. We also thank Andrew Lyne for providing the radio 
ephemeris and comments on the manuscript. Many useful discussions with Ed Fenimore 
and James Theiler are gratefully acknowledged. We are appreciative of the referee, F. 
Seward, for his many helpful comments in improving this paper. This work was performed 
under the auspices of the US Department of Energy and was supported in part by the 
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Fig. 1. — Epoch-folded lightcurves of PSR B1951+32 in three different energy bands. 
Two periods are shown for clarity. Phase is referred to the radio peak. In addition 
to instrumental and cosmic X-ray background, the DC component is mostly likely the 
contribution from the ROSAT point source and its associated X-ray nebula. 

Fig. 2.— The spectrum of PSR B1951+32 from 0.1 keV up to 10 GeV. The ROSAT data is 
for the X-ray point source including pulsed and unpulsed components. The dashed line is the 
solid line times 0.35, the reported pulsed fraction of this point source. The dotted line is a 
power law with a photon spectral index —3/2 passing through the COMPTEL data point in 
0.75-3.0 MeV. Upper limits are at a 2-a level for XTE/PCA, OSSE, and COMPTEL, and 3-ct 
for BASTE. References: EGRET - Fierro (1995); COMPTEL - Kuiper at al. (1996); OSSE 
- Schroeder et al. (1995); BATSE - Wilson et al. (1992); ROSAT - Safi-Harb, Ogelman & 
Finley (1995). 
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Table 1. Radio ephemeris of PSR B1951+32^ 



Validity interval 
Epoch, to 

^2000 
^2000 
V 

ii 

V 



(MJD) 
(MJD) 



(Hz) 

(Hz/s) 
(Hz/s/s) 



50057 - 50207 

50132.000000201 

19^52"^58^276 

32°52'40".68 

25.2963865292632 

-3.74022 X 10~^2 

2.81 X 10-22 



^provided by Andrew Lyne (1996, private communication) 
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